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Introduction: The spectral method for distin- 

guishing crustal from core-source magnetic fields has 
been re-examined, modified and applied to both a 
comprehensive geomagnetic field model and an alti- 
tude normalized magnetic map of Mars [1], These 
observational spectra are fairly fitted by theoretical 
forms expected from certain elementary classes of 
magnetic sources [2]. For Earth we found fields from 
a core of radius 3512 +/- 64 km, in accord with the 
3480 km seismologic radius, and a crust represented 
by a shell of random dipolar sources at radius 6367 +/- 
14 km, just beneath the 6371.0 km mean radius. For 
Mars we found only a field from a crust represented in 
same way, but 46 +/- 10 km below the planetary mean 
radius of 3389.5 km, and with sources about 9.6 +/- 
3.2 times stronger than Earth's [1]. 

It is remarkable that the same simple theoretical 
form should fairly fit crustal magnetic spectra for both 
worlds and return crustal-source depth estimates of 
plausible magnitude. Evidently, the idea of an ensem- 
ble of compact, quasi-independent, magnetized re- 
gions within these planetary crusts has some merit. 
Yet such estimates, at best a kind of average, depend 
upon both the observational spectrum fitted and the 
physical basis of the theoretical spectrum. 

Observational Spectra: The magnetic spectrum 
of a planet is the mean square magnetic induction 
configured in spherical harmonics of degree n, aver- 
aged over a centered sphere of radius r that contains 
the sources [3], 

n 

K(r) = (n + 1 )(a/r) 2n+4 2 [(g n m ) 2 + (/t n m ) 2 ]. (1) 
m = 0 

Here a is the radius of a reference sphere and (g n m , 
h a m ) are the Gauss coefficients of degree n and order m 
in a Schmidt-normal spherical harmonic expansion of 
the scalar potential V: B = -W. Observational spectra 
can be calculated from coefficients obtained via har- 
monic analyses of either measured data, or other mod- 
els fitted to such data. These spectra depend on the 
data fitted and modeling technique. 

For Earth, the magnetic held model we used has 
since been revised [4], updated with high precision 
0rsted data, and upgraded with high resolution 
Champ data [5]; at high degrees, the refined spectra 
appear softer than before. For Mars, there are now 
many fine models of the outstanding global MGS 
MAG/ER data [6 - 13]; these agree well through in- 


termediate degrees, but differ at hard to resolve high 
degrees. 

Theoretical Spectra: Consider a thin crust with 
compact, effectively dipolar, sources. If we expect 
dipole positions to be uncorrelated, random samples of 
a uniform distribution on a spherical shell, and vector 
dipole moments to be uncorrelated, random samples of 
a zero mean distribution, then our expectation spec- 
trum from an ensemble of such random dipoles on a 
shell of radius r x is [1] 

{ R n (a ) } ss = A x n (n + V 2 ) (n+l) (r x /a) 2n " 2 . (2) 

In spectrum (2), amplitude A x is proportional to 
the mean square moment of the sources and a cubic 
polynomial of n modulates the exponential attenuation 
function of rja < 1 . This spectrum increases with n at 
low degrees, peaks near degree 3/2/n(a/r x ), and falls of 
exponentially at high degrees. Curiously, the mathe- 
matical derivation of (2) [2] is analogous to that for 
magnetic change induced by compact hydromagnetic 
eddies near the top of a fluid conductor [14]. 

More realistic theoretical spectra, which allow for 
crustal thickness, oblateness and magnetization by a 
planet centered dipole, have already been derived and 
discussed, as have important spectral effects of later- 
ally correlated sources [1]. The latter were described 
via an ensemble of vertically and uniformly magnet- 
ized spherical caps; yet interesting complications arise 
from trial distribution functions for the sizes and 
magnetizations of extended sources. These are ad- 
dressed in hopes of improving source-depth estimates, 
resolving joint distribution parameters, and determin- 
ing whether or not the different crustal magnetic di- 
chotomies of the two planets, evident in the phase 
information, have distinct spectral signatures. 

Correlated Sources: The main effect of laterally 
correlated sources is to soften the spectrum at high 
degrees [1]. We tend to over-estimate source shell 
depth when this is omitted. To include this effect sim- 
ply, size and magnetization distribution functions for 
extended sources are recast as the characteristic half- 
angle v [/ 0 (hence diameter) and mean square total mo- 
ment { 7' 2 } for an ensemble of vertically and uniformly 
magnetized spherical caps on the shell of radius r x . 
The resulting theoretical spectrum is 

|R n (a)} = A (nil) [Z n ( ¥o )] 2 ( rja) 2a ~\ 


(3) 
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Amplitude A is proportional to { 7’ 2 } and, in terms 
of the Schmidt-normal associated Legendre polynomi- 
als _P n m (cos\|/) , Z n (\|/) = sin\|/ L ) n 1 (cosv|/)/[l-cost|/]. 

The physics of spectrum (3) is perhaps most easily 
seen via its finite Taylor expansion in e = 1 - cos\|/ 0 , a 
quantity proportional to the characteristic area of 
small source regions. To first order in small s, 

{R n (a)} « An 2 (n+V){rJd) 2 "- 2 

x [1 - (Vo 2 /4)7!(/!+1)]. (4) 

For small caps, and at moderate degrees n\\i 0 « 1, 
the partial derivatives of the logarithm of theoretical 
spectrum (4) with respect to amplitude A, shell radius 
r x , and cap half-angle \\i Q are approximately propor- 
tional to 1, n, and -n 2 , respectively. Separation of 
characteristic source diameter from amplitude and 
depth should thus be straightforward, unlike separa- 
tion of amplitude from layer thickness. The negative 
sign of the partial w.r.t. \j/ 0 describes the softening of 
the spectrum due to the small, but non-zero, area of 
the still unresolved sources. 

Indeed, with n\\i 0 « 1, ln[ 1 - (v|/ 0 2 /4)n(«+l )] is 
about -(i|/ 0 2 /4)«(h+ I ) and the logarithm of demodu- 
lated spectrum (4) is about 

ln[{R n (a)}/n 2 (n+l)\ ~ InA + 2(n - 1 )ln(a/r x ) 

- ( Vo 2 /4)«(«+l). (5) 

This is a parabola in n. Of course, spectrum (3) is 
positive, so approximate spectrum (4) and relation (5) 
require \\i a 2 « 4/n(n +\ )]. For high degrees and/or 
broader sources, higher order terms must be retained 
in the expansion of spectrum (3). For example, the 
second order correction to spectrum (4) is less than the 
first order term only for \| i 0 2 < (9.6) 1/2 /«(h+ 1). For 
degrees n < 90, this is so for source regions up to a 1.1 
degree half angle, or 130 km diameter, on Mars. 

Summary: Perhaps the simplest statistical model 
for compact sources in the magnetic crust of a planet 
is that of random dipoles on a shell. By carefully fit- 
ting a theoretical spectrum like (2) to an observational 
magnetic spectrum, we can estimate the characteristic 
amplitude and depth of such sources [1]. 

A more realistic model is needed to improve such 
estimates.. Perhaps the most important improvement 
is to include the small but non-zero areas of the source 
geologic structures. To this end. theoretical spectra 
have been developed for a simple class of laterally 
correlated sources. 

The results, here summarized simply via approxi- 
mate spectrum (4) and relation (5), yield spectra softer 


than (2). Moreover, comparison with previous spec- 
tral forms indicates that, in addition to source ampli- 
tude and depth, observational magnetic spectra from a 
planetary crust are likely telling us far more about the 
typical breadth of source regions than about typical 
source layer thickness. 

Estimation of three fundamental properties of a 
thin magnetized crust, the characteristic magnetiza- 
tion, depth, and breadth of its presumably small con- 
stituent source regions, boils down to a parabolic fit to 
the log of its cubically demodulated magnetic spec- 
trum. Moreover, the restrictive suppositions of small 
source regions and moderate harmonic degrees can be 
eliminated, though the computations, notably the for- 
mulas for the partials, become more complicated. 
Preliminary results from applications to observational 
spectra are to be presented and discussed. 
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